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ABSTRACT 

Massive high-redshift galaxies form in over-dense regions where the probabiUty of 
forming other galaxies is also strongly enhanced. Given an observed flux of a galaxy, the 
inferred mass of its host halo tends to be larger as its inferred redshift increases. As the 
mass and redshift of a galaxy halo increase, the expected clustering of other galaxies 
around it gets stronger. It is therefore possible to verify the high-redshift identity 
of a galaxy (prior to an unambiguous spectral identification) from the clustering of 
ot her galaxies arou n d it. W e illustrate this method for the massive galaxy suggested 
by iMobasher et al.l (120051 ) to be at redshift 2 ~ 6.5. If this galaxy were to exist at 
z '^ 6.5, there should have been a mean of ^ 10 galaxies larger than a hundredth of 
its mass and having z-band magnitudes less than ~ 25 detected as i-dropouts in the 
HUDF. We calculate an approximate probability distribution for neighbor galaxies and 
determine that there is less than a ^ 0.3% chance of detecting no massive neighbor 
galaxies. The lack o f other massive z ^ 6.5 galaxies in the HUDF image argues that 
the Mobasher et al.l p005 ) galaxy is instead a low redshift interloper. We generalize 
our results to other galaxy masses and redshifts. 
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1 INTRODUCTION 

The purported detection of a very massive galaxy, H UDF- 
JD2, at a redshift of 2~6.5 bv lMobasher et al.l l|2005h pro- 
vides a critical test for the standard paradigm of galaxy 
formation in a ACDM cosmology. However, the lack of an 
unambiguous spectral identification allows this galaxy to be 
a z ~ 2 interlo per in which obscuration by dust mim ics the 
Lyman break (jPunlop et al.ll2007l : [Charv et al.ll2007^ ■ If the 
galaxy is at z~6.5, the estimated mass of the halo contain- 
ing this galaxy is Mjd2 ~ 2 x 10^'^ Mq for a reasonable star 
formation efficiency, /^, of 10%. In the concordance cosmo- 
logical model, the expected number of such galaxies in the 
Hubb le Ultra Deep Field (H UDF) field-of-view is less than 
10~^ (jBarkana fc Loebl 120061 ). Since the expected probabil- 
ity of finding such a galaxy in the HUDF is extremely low, 
its redshift identification is of great importance for testing 
the standard cosmological model. 

Here we point out that, despite the low average abun- 
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dance of galaxies like JD2 (due to its large mass and high 
redshift), such a galaxy cannot exist alone. The large-scale 
over-density implied by its existence naturally results in 
neighboring halos over and above what would be expected 
from random fluctuations in the average galaxy population. 
The surrounding over-dense region behaves as if it is part 
of a closed universe, in which the formation of all galaxies 
occurs earlier. We approach this problem analytically in this 
Paper. In addition to the deeper fundamental understanding 
gained from such a treatment, the extreme rarity of objects 
like JD2 make a statistical analysis using numerical simula- 
tions difficult. In Sj2l we show how the excursion set formal- 
ism can be used to calculate an approximate probability dis- 
tribution for the number of neighbors around massive galax- 
ies. We then calculate, in par ticular, the expected cl ustering 
of bright galaxies around the IMobasher et al.l (|2005l ) galaxy 
in 321 and explore the dependence of our results on the star 
formation efficiency, duty cycle, and power-spectrum nor- 
malization in SJH In JlSl we generalize these results to other 
halo masses and redshifts. Finally, JS] summarizes our main 
conclusions. 

Unless otherwise noted, we assume a flat, ACDM model 
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for the universe wit h the WMAP3 cosmological parameters 
(|Spergel et al.ll2007l ). 



2 METHOD 

We assume the simpl e model for Lym an-break Galaxies 
(LBGs) considered by IStark et al.1 l)2007h . which associates 
LBGs with merger-activated star formation in dark-matter 
halos and includes suppression of the star formation ef- 
ficiency in low-mass halos by supernova feedback. In the 
model, the star formation duty cycle, enc, gives the frac- 
tion of halos that contain active star formation. This frac- 
tion has recently been calibrated by the measured luminosity 
function of LBGs at z ~ 6 to a best-fit va lue with 1 — a er- 
rors of £DC = 0.25to;o9 (|Starket al.ll2007h . Here we adopt a 
conservative value of eoc = 0.14, in accordance with our as- 
sumed star-formation efficiency of /* — 10% (which matches 
the fraction of Q,i, in stars today). The remaining fraction, 
1 — eoc, of halos at 2~6 will not be detected as i-dropout 
LBGs and may include a population of post-starburst galax- 
ies similar to JD2 itself. For our calculations, we assume no 
variation of the duty cycle over the redshift range of the 
selection function of the HUDF. 

Halos with active star formatio n do not constitute a fair 
sampl e of the total halo population. IScannapieco fc Thackerl 
(120031 ') show that these halos have undergone substan- 
tial accretion in their recent past giving them an ex- 
tra "temporal" bias. W hile the numerical simulations by 
IScannapieco fc Thackeij were done at z = 3, there is, as yet, 
no analytical method to predict this extra bias at higher 
redshift. Thus, our calculations for the number of neighbor- 
ing LBGs around massive galaxies at high redshift are lower 
limits that could be modified in the future with a better 
understanding of the evolution of the "temporal" bias with 

mass and redshift. 

According to the excursion set prescription (jZentnerJ 

120071 ). if the linear density fluctuations in the universe are 
extrapolated to their values today and smoothed on a co- 
moving scale R, a point whose over-density exceeds a critical 
value of Sc(z) « 1.686 D{z = 0)/D(z), where D(z) is the lin- 
ear growth factor at redshift z, belongs to a collapsed object 
with a mass M = (4/3) tt pcrtt R'^ if R is the largest scale for 
which the criterion is met, where pcrit is the critical density 
of the universe today. The critical value of the over-density, 
extrapolated to today from z = 6.5, is 5c{z — 6.5) ~9. 6. For 
a Gaussian random field of initial density perturbations, as 
indicated by WMAP3 measurem ents of cosmic microwave 
background (jSpergel et al.l 120071 ). the probability distribu- 
tion of the extrapolated and smoothed over-density, 5u, is 
also a Gaussian: 



MSR,SiR))d5-- 



exp 



^2nS{R) "'"" V 2S'(i?) 
with zero mean and a variance given by: 

dk 

2^2 



S{R) 



■k^ P(k), 



dS, (1) 
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where P (k) is the linear power-spectrum of density fluctu- 
ations today as a function of wave-number k, and kmax = 
1/R. Since equation ([2} is a monotonically decreasing func- 
tion of R (or Af ) , the smoothing scale can be uniquely spec- 
ifled by the variance of the over-density field smoothed on 



that scale. The critical threshold for collapse introduces a 
small correction to the probability distribution such that 
the distribution of Sr becomes: 

Q{5r, S{R)) = Qo{5r, S{R)) - Qo{2 5, - 5r, S{R)). (3) 

The conditional probability distribution of 5i on a scale 
specified by 5*1 given a value of ^2 on a scale larger scale 
specified by 5*2 < 5*1 is: 

Q{5i,Si\52,S2) = Q{Si - 52, Si ~ 82). (4) 

Using Bayes Theroem, the conditional probability of 82 on 
a scale 52 given Si on a smaller scale specified by Si > S2 
is: 

Q{52,S2\5i,Si)dSi(xQ{5i,Si\52,S2)Q{52,S2)d52 (5) 

with a constant of proportionality such that the integral of 
equation ([5} is unity. Setting 5i — 5c{z) and 5*1 — S(Afi), 
where Afi is the host halo mass of a detected massive galaxy, 
gives the probability distribution of ^2 = 5r on any comov- 
ing scale R > Ri, due to the presence of that galaxy, where 
_Ri is the radius corresponding to Afi. 

However, the excursion set formalism calculates the col- 
lapse of objects (a nonlinear effect) by considering the be- 
havior of the linear over-density field extrapolated to the 
present day. This method functions entirely in Lagrangian 
coordinates (which move with the flow) and does not take 
into account how the over-density field changes in the quasi- 
linear regime. Obviously, when a region collapses, matter is 
pulled in from the surrounding region to fill the void. Thus, 
if we assume the existence of JD2 at z = 6.5, the material 
in the rest of the HUDF at that redshift would have started 
outside the region earlier in the universe's history. 

We denote the Lagrangian radius of a region as Rl- 
Early in the history of the universe, before the region begins 
to collapse, this radius is equal to the radius of the region in 
Eulerian coordinates (which do not move with the fiow) . As 
the region collapses, the Eulerian radius shrinks, while Rl 
remains unchanged. We denote the final Eulerian radius of 
the region at the redshift at which it is observed as Re- 

The extent of the collapse depends on the magnitude 
of the over-density in the presence of the massive galaxy 
whose probability distribution we have just calculated in 
Lagrangian coordinates. The more over-dense the region, 
the larger it would have to be initially to collapse to the 
same value of Re ■ Similarly, a lower value of the over-density 
would mean that the material inside Re came from a rel- 
atively smaller Lagrangian size. We would like a mapping, 
then, between the comoving Eulerian size of a viewed re- 
gion (such as the HUDF), and the comoving Lagrangian 
size of the region from where the same material originated 
in the early uni verse. This can be o btained via the spherical 
collapse model (|Mo fc Whitelll996l ). A spherically symmet- 
ric perturbation of Lagrangian radius Rl and over-density 
5l > collapses to a sphere of comoving Eulerian size Re 
at redshift z given by: 



3 1 - cos 6* D(z = 0) „ 
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(6) 



(7) 



For a fixed value of Re, there is a one-to-one relationship 
between Rl and the value of Sl that collapses Rl to Re- 
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Figure 1. The upper panel shows the probability distribution of 
the over-density 5 (extrapolated to z = 0) due to cosmic vari- 
ance (Eq. ISjl in a spherical region that collapsed to the size of 
the HUDF assuming the existence of a 2 X 10^'^ Mq halo con- 
taining JD2 at z = 6.5. The central panel shows, for each value 
of the over-density, the resulting number of LBGs in halos with 
mass above 2 X 10^^ Mq expected in this region. The final prob- 
ability distribution of the number of LBGs in the region above 
2 X 10^^ Afo is plotted in the bottom panel, taking into account 
both cosmic variance and Poisson fluctuations. 



In the presence of a massive galaxy, the probability dis- 
tribution of S (in the Lagrangian sphere that collapsed to 
Re) can be computed by considering the possible histories of 
the region Re having collapsed from difTerent possible Rl 's 
weighted by the probability of the corresponding value of 5^ 
in each Rl ■ These weights are given by equation (O . The re- 
sulting probability distribution of 5 (dropping the subscript 
L) "seen" in a fixed Re can be expressed generally as: 



dP{5\Mi) 
dS 



oc Q{5, Rl{S, Re)\Sc{z), R{Mi)), 



(8) 



where again the constant of proportionality is set so that 
J{dP{5\Mi)/d5)d5 = 1. 

Now that we know the distribution of over-densities in 
which the massive galaxy sits, we can easily calculate the 
expected number of neighbor galaxies seen in Re in the 
presence of each value for the over-densi ty, N(S\Mi), and 
thus, the probability distribution of N. iBarkana fc Loebl 
(I2OOJ) calculate the mass function of objects in a region 
of fixed over-density by combining the Sheth-Tormen and 
Press-Schechter prescriptions in the regimes for which each 
best fits results from numerical simulations. Given the rela- 
tionship between S and Rl for a given Re, the over-density 
determines the Lagrangian scale size. We find 



N{S\Mi) = {V(S) - V{Mi)) eoc 



dm 



{m,z,5) dm 
(9) 



dnuas , jx dnsT , ,fps{Sc{z)-S,S{m) — S{5)) 
-(m,z,d) — — (m,z)- 



dm 



dm 



}ps{5.{z),S{m)) 



(10) 

where dnsr/dm, is the Sheth-Tormen mass function, V{S) is 
the volume enclosed by the Lagrangian radius specified by 
5, and 



fpsiS,iz),S)dS- 



Sc{z) 



MSciz),s)ds 



(11) 



53/2 

is the mass fraction at z contained in halos with mass in the 
range corresponding to (S, S + dS) . The mass limit rrimax is 
the mass enclosed by a sphere of today's critical density with 
radius Rl corresponding to 5. We ignore the probability that 
Rl, while containing an over-density 5 < Sc{z), might be 
part of a larger collapsed region with 5 > Sc{z). This is a 
good assumption given the unlikely occurrence of 5 = Sc{z) 
on the scale of Mi in the first place. 

The resulting probability distribution of iV, due to cos- 
mic variance, is given by: 

P{N\Mi) _ dP{S{N)\PIi) dS 

dN' 



dN 



dS 



(12) 



where 5{N) is the inverse of equation Q and dS/dN is its 
derivative. Poisson fiuctuations contribute additional vari- 
ation in the actual number, A*', of neighbor galaxies. The 
probability, P{N\Mi), of each discrete value of A'^ in the 
presence of a mass Mi galaxy at high redshift is obtained 
by convolving a discretized version of equation (|12|) with the 
Poisson distribution in the following way: 
00 
P(iV|Afi) = Y. P{N\Ml)Ppo^sso„{N,N), (13) 

JV=0 

where 



son\f^, Aj 



P{N = 0|Mi) 



and for N > 0, 
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P{N'\Mi) 
dN' 



dN, 



dN'. 



(14) 



Equation l|13|l can be compared to galaxy counts in surveys. 
In particular, if no galaxies in halos above some minimum 
mass, rrimin, are seen as neighbors to another galaxy in a 
halo of mass Mi, then the quantity 1 — P{N = 0|Mi) is 
approximately the confidence by which we can rule out ei- 
ther the existence of a halo of mass Afi or the cosmological 
model. 



3 HUDF-JD2 

We show results for neighbor i-dropouts around JD2 in the 
HUDF. If JD2 is indeed a massive galaxy at z ~ 6.5, then 
there should be many massive galaxies visible around it. Fig- 
ure[T]shows the distribution of 5 given by equation ((Sjl inside 
the Lagrangian patch that collapses to an angular Eulerian 
size of '^IIS", enclosing an area on the sky roughly equiv- 
alent to the HUDF field-of-view. This angular scale corre- 
sponds to a comoving size of ~4.7Mpc at z — 6.5, which 
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Figure 2. The dependence of the number of neighbor galaxies 
on the star formation efficiency and the duty cycle assuming a 
fixed value of 2 X 10^^ Mq for the stellar mass of JD2. The upper 
two panels show the mean number of neighbor LBGs at z = 6.5 
with a host halo mass above 0.01 X Mju2 = 2 X 10^^ Mq in 
the HUDF due to the presence of JD2 and the probability of 
detecting none of these objects as a function of the duty cycle. 
The assumed one-to-one relationship between the star formation 
efficiency and the duty cycle given in equation I I15I I is plotted in 
the third panel, while the lower panel shows how the host halo 
mass of JD2 depends on the duty cycle through its relationship to 
the star formation efficiency. The vertic al dashed line deno tes the 
best-fit value of the duty cycle given bv lStark et al.l 1 120071 ). while 
the dotted lines indicate the 1 — ct bounds. Since the validity of 
equation I I15I I cannot be verified beyond the range of the 2 — a 
contour, we truncate the plot near its lower boundary at euc of 
~0.1. 



6.5, encloses 
0. Also shown 



assuming an extrapolated over-density of 5 
within a spherical radius a mass of ~ 10^"* M^ 
are the expected number of neighbors inside that radius, 
given each value of the over-density, and the resulting prob- 
ability distribution of the number of neighbors that includes 
both cosmic variance and Poisson fluctuations. 

While the HUDF should be sensitive to LBGs in hosts 
as small as ~ 2 x 10^" Mq , the number of random galax- 
ies expected in such hosts without correlations from JD2 is 
comparable to the number of excess neighbors that result 
from these correlations. This creates some ambiguity in de- 
tecting the excess over the background. On the other hand, 
while the mean number of uncorrelated LBGs in halos as 
large as ~ 2 x lO'^^ Mq is orders of magnitude smaller than 
the excess due to JD2, the probability of detecting no such 
neighbor is not small enough for a lack of a detection to be 
meaningful. 

Thus, we consider neighbor LBGs in halos with masses 
above 2 x lO^'^ Mq. For a duty-cycle eoc = 0.14 and star 
formation efficiency /^ = 0.1, these galaxies have lumi- 
nosities at 1500 A above ~2 x lO^^ergss"^ Hz~^ or z-band 
magnitudes at z ~ 6.5 less than ~ 25. The mean abun- 



Figure 3. The dependence of the number of neighbor galaxies on 
the normalization of the matter power spectrum ag. The upper 
panel shows the mean number of neighbor LBGs at 2 = 6.5 with 
a host halo mass above 0.01 x Mj/32 in the HUDF due to the 
presence of JD2, while the lower panel indicates the probability 
of detecting no such objects. The vertical dashed line denotes the 
WMAP value used in the rest of the Paper. 



dance of such galaxies in the absence of JD2 is much less 
than unity, an d indeed, no such obj e cts ha ve been detected 
in the HUDF (|Bouwens et al.ll2006l . [2007h . Figure [T] shows 
that there should be a mean of ~ 10 very bright i-dropout 
LBGs in hosts with masses larger than within an angular 
Lagrangian radius of '^ 175" of JD2, where this is a lower 
limit due t o the fact that LBGs should be more clustered 
than halos (|Scannapieco fc Thackeiil2003l ). The probability, 
given by equation H13p , of detecting no such galaxies in this 
region is P{N = 0[Mjd2)~3 x lO""*. Thus, we can rule out 
JD2 at redshift z ~ 6.5 with 99.7% confidence. Integrating 
the distribution, we find a less than a 5% chance of detecting 
fewer than 3 very bright neighbor galaxies. Either JD2 is at 
z^2, as allowed by spectral fits, or there is a problem with 
our assumed cosmological model. 



4 PARAMETER DEPENDENCE 

So far, we have assumed a star formation efficiency of /* = 
10%, a duty cycle of enc ~ 0.14, and WMAP3 cosmological 
parameters with as = 0.776. We now explore how our results 
depend on these parameters. 



4.1 enc and /* 

Given a stellar ma ss of 2 x 10^^ Mq for JD2 
IjBarkana fc Loebl 120061 ) . changes in the star formation 
efficiency will affect the assumed host halo mass. More 
efficient star formation will, therefore, result in fewer 
neighbor halos due to the lower corresponding halo mass. 
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Figure 4. Results for future detections of massive galaxies at 
fixed zi as a function of host halo mass Mi. The upper panel 
shows the mean number of neighbor LBGs with a host halo mass 
above 0.01 X Mi within the angular Eulerian distance correspond- 
ing to the HUDF, while the probability of finding no such LBGs 
in the HUDF is plotted in the lower panel. The dotted, solid, 
and long-dashed lines denote values for 21 = 5, 6.5, and 8, re- 
spectively. The vertical dashed line denotes the host halo mass of 
JD2. 




Figure 5. Results for future detections of massive galaxies as a 
function of redshift z. The upper panel shows the mean number 
of neighbor LBGs with a host halo mass above 0.01 X Afi in 
the HUDF, while the probability of finding no such LBGs in the 
HUDF is plotted in the lower panel. The solid, dotted, and long- 
dashed lines denote values for Mi = 2 X 10^^ Mq, 10^^ Mq, and 
2 X 10^^ Mq, respectively. The vertical short-dashed line denotes 
z = 6.5, the redshift of JD2. 



Meanwhile, changes in the duty cycle will affect the fraction 
of such halos that are seen as i-dropouts. Yet, these two 
parameters are no t constrained independently. Figure 4 of 
IStark et al.l ([20071) shows very narrow likelihood contours in 
log-parameter-space. As a rough approximation, then, we 
assume that the range allowed by the luminosity function 
fitting spans a power-law relationship between the star 
formation efficiency and the duty cycle. Using the best-fit 
parameter values and the extremes of the 1 — a contour, 
we fit the dependence with a least-squares regression (in 
log-space) and find a relationship given by: 



f*^Ae[ 



DC 



(15) 



where A = 0.264 and f3 = 0.378. 

The mean number of neighbor i-dropouts to JD2 in the 
HUDF with host halo masses greater than 0.01 AIi, denoted 
(A'^), is plotted in Figure[2]as a function of eoc- Also shown 
is the effect on the probability of detecting no such galaxies, 
given by setting A^ = in equation (|13p . As shown, the effect 
of varying parameters on (A'^) is only modest. When consid- 
ering the range of 1 — ct errors on eoc, {N) varies by only 
a factor of ~ 2. The correlation between the star formation 
efficiency and the duty cycle given by equation (|15|l . causes 
these parameters to moderate each other's effect on the ex- 
pected number of neighbors. While the value of P{N = 0) 
varies more significantly, our choice of parameters through- 
out the rest of the paper gives conservative values for both 
{N} and P{N = 0). The best-fit value of eoc = 0.25 derived 
by Stark et al. (2007) results in an even smaller probability 
of detecting no bright neighbor dropouts if JD2 is at 2 = 6.5. 



4.2 erg 

The high over-density in the region surrounding a massive 
galaxy at high redshift results from the large density fluctu- 
ation required to produce such a galaxy. If the over-density 
must reach 6c{z) on the scale corresponding to the size of 
the galaxy, then on a smoothing scale only a little larger, the 
over-density could not have been very low, since the contri- 
bution from the intervening scales is a Gaussian about zero 
with a standard deviation much less than Sc{z). However, 
varying the normalization of the matter power-spectrum will 
change the standard deviation of this Gaussian contribution. 
Reducing the value of as will decrease the contribution from 
these intervening scales and cause the over-density around 
massive galaxies to be larger. This will result in more neigh- 
boring galaxies on average and a lower likelihood of detect- 
ing none. Conversely, increasing as will boost the possible 
contribution from intervening scales, shift the probability 
distribution of S in regions around massive, high-redshift 
galaxies toward lower values, and result in lower average 
number of neighbors and a greater chance of detecting none. 
We explore the dependence of (A*') and P{N = 0) on as 
quantitatively in Figure [S] The behavior is just as expected. 
The dependence is relatively strong given the wide range of 
proposed values for as- However, P{N = 0) reaches only 
~ 1% at as = 0.9. 
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5 POTENTIAL FUTURE DETECTIONS 

As new surveys conti nue to probe for ever larger galaxies at 
ever h igher redshifts IIMcLure et alJbood : iRodighiero et alj 
I2OO7I : IWiklind et al.1 l2007f ). it is important to generalize 
our results to potential future detections of other mas- 
sive, high-redshift galaxies in the future. We again assume 
[eDC, f*,o-s) = (0.14,0.1,0.776) and calculate the expected 
number of neighbors as for JD2 but now allowing the mass, 
Ml, and redshift, zi, of the halo hosting the detected galaxy 
to vary. 

Figure [3] shows the mean number of neighbor LBGs, 
(A'^), with host halo mass M2 > 0.01 x Mi that are within an 
angular Eulerian separation 6 = 115" from a central galaxy 
with a host halo mass of Mi situated at a fixed zi as a func- 
tion of Ml. While the uncorrelated, average abundance of 
lower mass halos is larger than that of higher mass halos, 
the plot demonstrates that the average number of neighbors 
begins to increase as a function of Mi, for zi = 6.5 and 8 
above Mi ~ 5 x 10^^ Mq due to the nonlinear increase in 
halo clustering with mass. Below this value, the increasing 
uncorrelated abundance of objects with lower mass causes 
the number of neighbors to increase with decreasing Mi . At 
zi = 5, (A'^) decreases as Mi increases in the entire range 
considered. The correlative effects of even a very massive 
halo are dominated by the uncorrelated behavior at this 
redshift. The mean number of neighbor LBGs within 115" 
is greater than unity for all values of Mi and zi considered. 

The probability of detecting no such neighbors, P{N = 
0), is also plotted in Figure U) Its dependence on Mi is as 
expected; P{N — 0) increases as {N) decreases. 

In FigureO {N) and P{N = 0) are plotted as functions 
of zi at fixed Mi and 6 = 115". While the plots of (A^) and 
P{N — 0) increase and decrease, respectively, with increas- 
ing 21 for Ml = 10" Mq and 2 x lO" Mq due to increased 
clustering with the detected halo, this is not indicated for 
Ml = 2 X 10^^ Mq. As shown in Figure H d {N) jdMi < 
at Ml = 2 X 10^^ Mq for all values of zi considered indi- 
cating that clustering is less important for these halos and 
their neighbors. This is also seen in the larger values of (A^) 
for Ml = 2 X 10^^ Mq compared to lO" M© . The increase 
is due to an uncorrelated population of halos with mass 
2 X 10^° Mq < M2 < 10" Mq. For each set of parameters, a 
higher value (A'^) corresponds to a lower value of P{N — 0), 
as expected. 



6 DISCUSSION 

We have presented a clear signature for the existence of a 
massive galaxy at high redshift. Due to large correlations 
among massive halos at high redshifts, once such an object 
is found, it is unlikely to be alone. On the contrary, the more 
massive the halo and the higher its redshift, the greater the 
number of neighbors it has. This signature is helpful in sit- 
uations where only photometric data is available or when 
the spectroscopic redshift identification is ambiguous. For 
a galaxy of a given observed flux, the inferred halo mass 
increases as its suggested redshift increases. It therefore be- 
comes easier to rule out a higher-redshift hypothesis, using 
dropout techniques to locate the neighboring LBGs that are 
expected at a similar redshift. 



We derived a probability distribution for the number of 
neighbors around a massive galaxy at high redshift, which 
led to calculations of the mean number of such neighbors 
and the probability of detecting no such halos. 

In the case of JD2, the number of neighbor galaxies ob- 
served in th is way can distin g uish b etween the z~6.5 inter- 
pret at ion_ofJMobashe£^tayjl2005|}_Md_tlw^ of a galaxy at 
z ~ 2 (|Dunlop et al.ll2007l : ICharv et al.ll2007^ . If at z = 6.5, 
the mean number of excess i-dropouts in the HUDF with z- 
band magnitudes less than ~ 25 is predicted to be ~ 10. The 
lack o f such bright i-dropouts in the HUDF (jBouwens et al.l 
|2006| . 120071 ') implies that JD2 cannot be at such high red- 
shift with 99.7% confidence. Thus, future detections of mas- 
sive galaxies at high redshift could be verified by looking for 
bright, massive neighbors. 

Additional research into exactly how LBGs form in 
dark matter halos and how biased they are compared to 
these halos would firm-up the predictions for the probabil- 
ity distribution of neighbors. Such work would undoubt- 
edly include various feedback processes and environmen- 
tal effects on LBG formation. Previous study indicates 
that LBGs are more cluster ed than dark matter halos 
flScannapieco fc Thackeiil2003l ) , and our results for the num- 
ber of neighbors are, thus, lower limits on the true values. 
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